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We present a field theory approach to study changes in local temperature due to an 
applied electric field (the electrocaloric effect) in electrolyte solutions. Steric effects 
and a field-dependent dielectric function are found to be of paramount importance 
for accurate estimations of the electrocaloric effect. Interestingly, electrolyte solu¬ 
tions are found to exhibit negative electrocaloric effects. Overall, our results point 
toward using fluids near room temperature with low heat capacity and high salt 
concentration for enhanced electrocalorics. 


There has been a renewed interest in developing caloric materials 


technologies [l|, Q] for various refrigeration applications. The caloric materials undergo 
reversible thermal changes under the influence of an applied field, which can be mag¬ 
netic, electric or mechanical in nature. These thermal changes due to magnetic field, 
electric field and mechanical stresses are known as the magnetocaloric, electrocaloric and 
mechanocaloric effects, respectively. Thermodynamic description of these changes was pro¬ 
vided by Thomson}^ and the changes are results of variations in entropy of the system 
under the influence of an applied field. The magnetocaloric effect is already used to reach 
temperatures in the milliKelvin (mK) range and is in the stage of being commercialized for 
household refrigeration. In contrast, search for novel materials that can achieve the so-called 
colossal or giant electrocaloric effect is currently a topic of extensive research. Historically, 
ferroelectric materials[l-4, 6|, which are crystals with net polarization in the absence of any 
external applied electric field, have been studied extensively for the electrocaloric effect and 
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have shown thermal changes as low as 0.003 K near room temperatnre and as high as 31 
K based on the operating temperature and applied electric heldj^. It is to be noted that 
most of the materials studied for the electrocaloric effects are in the solid state except some 


polymeric films 


[ll-y,[^, which are considered viscoelastic, and liquid crystalline fluids js] that 
have shown giant electrocaloric effects in thin him geometries. 

Thermodynamic description of the electrocaloric effect in an adiabatic system relies on 
the fact that changes in entropy resulting from application of an electric held must be zero. 
As entropy can be modihed by varying either temperature (T) or diherence in the surface 
potentials of the electrodes (lA), we consider inhnitesimal changes in entropy {AS) for a 
system at an initial temperature of T = Tq and the potential diherence Vg = Vq undergoing 
inhnitesimal changes in the temperature (AT) and the potential diherence (Aldj) so that 

dS 


AS = 


dS 

&T 


T=To, 

Vs=Vo 


AT+ 


dVg 


T=To, 

Vs=Vo 


AVg. 


( 1 ) 


For adiabatic changes, AS = 0 and noting that T ^ = c„(T = Tq, Id, = Vq), 

where is the volume heat capacity and depends on the initial temperature and the potential 
diherence, we can write 


'AT' 


e/ks 

'dS' 


T=To, 

Vs=Vo 

c,{T = To,Vg = Vo) 

d'ipg 


T=To, 

, _ eVp 


( 2 ) 


Here, we have dehned 'ijjg = eVg/k^TQ, e is the charge of an electron and ks is the Boltzmann 
constant so that e/ks = 1.16 x 10'^ K/V. It is to be noted that electrostrictionjs] ehects 
leading to changes in the volume of the liquids are not taken into account here and form the 
basis of multi-caloric materials exhibiting electrocaloric and mechanocaloric/elastocaloric 
ehects. This is an interesting direction for future research. Eq. [2] provides three insights. 
First, it is clear that the changes in temperature resulting from changes in the potential 
diherence are inversely proportional to the volume heat capacity of the material. Hence, 
huids with low heat capacity are preferable candidates for enhanced electrocalorics. Second, 


insight is obtained from the use of thermodynamic rules stating 


10[ | that entropy must 


in¬ 


crease with an increase the temperature i.e., > 0. This implies that the dimensionless 

quantity eld/fc^To, ratio of the electrostatic energy of a unit charge to the thermal energy, 
is the relevant variable. In particular, sign of changes in the temperature (i.e., increase or 
decrease) with an increase in the surface potential depends on the changes in entropy with 
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respect to eVg/ksTo. Third, the length scale of the region undergoing changes in temperature 
is determined by the volume undergoing entropic changes. 

Larger entropic changes resulting from small changes in the potential difference are re¬ 


quired for enhancing t 
are expected in liquids 


lie e 


ectrocaloric effect (cf. Eq. [2]). As larger entropic changes 


71, Isl] than solids in the presence of an external field, we have fo¬ 
cused on a theoretical description of the electrocaloric effect in electrolyte solutions. We 
use Eq. |2] and entropic changes computed using field theoryjh, 6| to study the elec¬ 


trocaloric effect in planar double layer systems 


layer can be constructed 


MM 


effects due to dielectric saturation 
of ions 0, y, 12^ 


19 

Q1 




The free energy of the double 


2C]with different approximations including various 
^ll, l22| , finite polarizability of ions 24j-l26| , finite size 
. y |. ion adsorption-desorption equilibrium M] and image charges jl4. 321. 


This allows systematic investigations into roles played by different factors in affecting the 
electrocaloric effect and pave the way for rational design of enhanced electrocaloric fluids. 
Another motivation in studying such a system lies in the need for an improved theory for the 
electrolyte solutions in strong external fields, where crowding and dielectric saturation effects 
are important and a larger electrocaloric effect is observed for viscoelastic materials such as 


polymer 0 films and liquid-crystalline solutions|8|]. Furthermore, novel technologies 33|, [34 1 




for extracting energy by mixing fresh river water with saline ocean water can benefit from 
an improved theory for the electric double layer. These technologies are based on the well- 
known fact that an electric double-layer acts as a capacitor and salt concentration plays a 
key role in dictating its capacitance. Operating t emp erature has been shown to play a key 
role in affecting the energy that can be harvested 3J] using these technologies. 

We use a microscopic field theory approach to study planar double layer systems (see the 
Supporting Information). In particular, we consider two parallel plates having surface area 
A each, separated by distance L and immersed in an electrolyte solution containing equal 
number density (= pc,b) of positive and negative ions along with ps^b as the number density of 
solvent molecules. The plates are assumed to have uniform surface charge densities (number 
of charges per unit area), ai and ct 2 and the corresponding surface potentials are Vo,i and 
Vo ,2 (in units of Volts), respectively. Surface potentials and charge densities are related to 
each other by electrostatic boundary conditions and depend on the mechanisms by which 
the plates acquire the surface charge. These relations can be formally derived by consider¬ 
ing different mechanisms for charging. We take molecular volumes of the solvent, positive 
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and negative ions to be Vs,v+ and n_, respectively. Noting that theoretical description of 
polarization nnder an external electric held and strong electric helds are pre-reqnisites for 
developing theory for the electrocaloric effect, held dependent dielectric and steric ehects 
resnlting from hnite sizes of ions and solvent molecnles are inclnded in onr model. In this 
work, we have bnilt a minimal model that can captnre the nnderlying physics based on treat¬ 
ing each solvent molecnle as an electric dipole of length pg occnpying molecnlar volnme Vg. 
Finite polarizability of ions and solvent molecnles are not considered in this work. However, 
the cnrrent formalism can be extended to take into acconnt the ehects of polarizability. We 
have nsed the theory to stndy the electocaloric ehects in non-overlapping donble layers (i.e., 
single donble layer systems) so that Vo,i = Vo and Vo ^2 = 0 i.e., conditions of constant snrface 
potentials are considered in this work so that the potential diherence Vg = Vq. Parameters 
are chosen to describe water molecules (such as the dipole moment pg). Furthermore, in 
these model calculations, we have considered symmetric ions and solvent molecules so that 
Vg = V-^- = V- = and ignored the asymmetry in sizes of the molecules. The size parameter 
a is chosen so that the density of pure water is reproduced i.e., l/a? = 1 gm/cm^. 

Typical free energy changes (AF) of the double layers (with respect to the electrolyte 
solution in the absence of applied surface potential ) are shown in Figure [11(a) for diher- 
ent values of eVo/^s^o cind temperature (Tq) ranging from room temperature to near the 
boiling point of water. The free energy changes are negative for the entire parameter range, 
which is in qualitative agreement with the predictions of the standard Poisson-Boltzmann 
(PB) approach (i.e., ignoring held-dependent dielectric and steric effects) and the modi- 
hed Poisson-Boltzmann (MPB) approach (i.e., ignoring held-dependent dielectric effects) 
(cf. Eqs. 47 and 43, respectively, in the Supporting Information). Also, larger free energy 
changes are found with an increase in the temperature due to increased entropic contri¬ 
butions shown in Figure [11(b). Furthermore, an increase in the free energy changes with 
an increase in the surface potential is also in qualitative agreement with the PB and MPB 
approaches. Corresponding entropic changes (AS* = —{dAF/dT)n,fl = AL being the total 
volume) such as those shown in Figure [11(b) dictate the electrocaloric effect. 

As the free energy and entropy changes per unit area are computed using the held theory, 
we rewrite Eq. [2lto calculate the electrocaloric ehect so that 


AT 


e/kB IdAS/Aks 


( 3 ) 


AVg T=To 


J- —J-K)-, 
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where A = A/a? and = c^/Aks is the rescaled heat capacity of the electrolyte solution 
in the presence of applied electric held. Formally, it can be written as = c^iT = Tq, Vs = 
0) + c„(T = To, Vs = Vo) so that c„(T = Tq, Id; = 0) = T 


dSh/AkB 

dT 


T=To, 

1ps=0 


is the rescaled 


heat capacity of the reference homogeneous electrolyte solution having as its entropy 


and c^(T = Tq, 14 = Vo) = T 


dAS/Akg 

dT 


T=To, 


accounts for additional contributions due 


1ps = 


eVp 




quantity 


to the applied electric held. It is to be noted that in Figure [Hb), surface potentials 
and temperature are varied simultaneously due to the variation of 'ipo = eVo/ksTo and the 
aAsjAkg extracted from Figure [U^b) using the formal relation = 

— ^ dAS/AkB ^ calculating the electrocaloric ehect, we have taken 

p,=Po 1 ™ -It=To 

c^(T = To, 14 = 0) = Q.QLdi/a corresponding to molar heat capacities of water to be 3.0i? 
(taken to be independent of temperature) and 3/2i? for each type of ion treated as an 
ideal gas 10| in the homogeneous phase, where Ldi is the thickness of the double layer and 
R = ksN^ is the universal gas constant so that is the Avogadro’s number. It is to be 
noted that Ldi naturally sets the length scale of the region undergoing changes in temperature 
as a result of the electrocaloric ehect. For the numerical estimates, we have dehned Ldi/a 
as the distance from the electrode after which counterion and coion densities approach their 
bulk values, pc,b, 

From isothermal changes in the entropy in Figure [11(b), it is clear that entropy of the 
double layer increases with an increase in the surface potential (i.e., (?AS'/54s > 0). Such 
an increase in the entropy is in qualitative agreement with the predictions based on PB 
and MPB approaches (see Eqs. 48 and 46,respectively, in the Supporting Information). As 
per Eq. [3l this should lead to a decrease of temperature with an increase in the surface 
potential i.e., a negative electrocaloric ehect is expected. Indeed, such a behavior is observed 
in Figure |2]for diherent initial temperatures and salt concentrations in the bulk. Figure [2] 
provides the magnitude of the electrocaloric ehect. As an example, consider an electrolyte 
solution containing 1.0 M monovalent salt with an electrode at surface potential of 0.2 V 
at Tq = 303 K (near room temperature). For this particular system, AT/A14 ~ —0.8 K/V 
is determined from Figure [2](b) so that a temperature decrease of 0.16 K is predicted. It 
is found that magnitude of AT/A14 is dependent on the initial temperature and the salt 
concentration in the bulk. In particular, the magnitude decreases with an increase in the 
temperature and increases with an increase in the salt concentration. The decrease in the 


















6 



(a) (b) 

FIG. 1: (a) Changes in the free energy and (b) entropy as a function of applied surface potential 
(Vb) and temperature for an electolyte solution containing 0.1 M monovalent salt. Legends show 
the values of eVo/^s^o 



FIG. 2: (a) Effects of initial temperature (Tq) and (b) the bulk salt concentration (so that Pc^ = 
0.6023cs (nm)“^ and Cg is in moles per litre (M)) on the electrocaloric effect. The left panel 
corresponds to = 0.1 M and the right panel corresponds to Tq = 303 K. Parameter for double 
layer thickness Ldi/a is found to be 3,4 and 6 for Cg = 1.0,0.5 and 0.1 M, respectively. 
























7 


magnitude with an increase in the initial temperature is a direct outcome of an increase in 
the heat capacity of the double layer with an increase in the applied surface potential, as 
evident from Figure [^b). The increase in the magnitude of the electrocaloric effect with an 
increase in the bulk salt concentration results from a decrease in thickness of the double layer 
{Ldi). Furthermore, larger free energy and entropic changes are found with an increase in the 
bulk salt concentration, as shown in Figure 1(a) in the Supporting Information. It should 
be noted that qualitatively the same effects are predicted by the PB and MPB approaches, 
where the free energy and entropy changes increase as ^Pc,b- However, quantitatively, the PB 
and MPB approaches digress from the numerical results due to errors made in predicting 
the free energy changes. To demonstrate this point, we have shown a comparison of the 
free energy changes for the same system, estimated using the PB, MPB and the numerical 
calculations in Figure ll](a). It is found that the PB approach is off by factors of 10 — 100 
for eVo/ksTQ > 10. In contrast, the MPB approach corrects for some of the errors made in 
the PB approach but it still deviates from the numerical results by a factor of 3. 



(a) 


(b) 


FIG. 3: (a) Comparison of the free energy changes computed using the PB approach (i.e., ignoring 
field-dependent dielectric and steric effects), the MPB approach (i.e., ignoring field-dependent 
dielectric effects) and the numerical calculations for Cg = 0.1 M, Tq = 303 K. (b) Computed surface 
charge density as a function of applied surface potential, estimated using the PB, the MPB and 
the numerical calculations, for different bulk salt concentrations at Tq = 303 K. The solid lines 
correponds to the analytical relation presented in the text. 




The relative importance of the field-dependent dielectric and steric effects in predicting 
structure of the double layer and resulting changes in the free energy can be assessed by 
comparing plots showing surface charge density as a function of the surface potential as 
predicted by the PB, MPB and numerical approaches (Figure 111(b)). The surface charge 
density in the PB approach is given hy ai = a = ^yj^^j^sinh where \Zc\ is the 

valency of ions (= 1 for monovalent ions), Ibo = /ATieoksT so that eo is the permittivity of 
vacuum, = 1 -|- 4:7ilBoPlPs,b/^ is the relative permittivity of the homogeneous electrolyte 
solution so that is the solvent density and k = {SttIboIZcI'^ pc,b /is the inverse 
Debye screening length. The PB and MPB approaches predict a monotonic increase in the 
surface charge density with an increase in the surface potential, as shown in Figure 111(b), 
without showing any sign of saturation, leading to unphysical surface charge densities. The 
numerical calculations show agreement with the PB and MPB approaches for eVo/ZcsTo < 
2 — 4 depending on the salt concentration and deviate strongly for higher surface potentials 
exhibiting saturation and a decrease in the surface charge density. This is an outcome of 
dielectric saturation leading to lowering of dielectric function near the surface, ignored in the 
PB and MPB approaches. The decrease in the surface charge density with an increase in the 
surface potential (i.e., negative differential capacitance[3a,|3^, where differential capacitance 
= da/dVo) hints at the breakdown of the one-dimensional uniform charge density model used 


here and plausible onset of in-plane charge density waves 


35j. 


In conclusion, we have presented a held theory approach for studying electrocaloric effects 
in planar double layer systems. Two key ingredients of the theory are the consideration of 
steric effects and dipolar interactions resulting from polar solvent molecules. Although 
the theory is general, in this work, we have presented calculations for aqueous solutions 
containing monovalent salt ions. It was shown that the electrocaloric effect in planar double 
layer systems is negative, i.e., the temperature of the double layer should decrease with an 
increase in the applied surface potential. The magnitude of the electrocaloric effect depends 
on the initial temperature of the solution and the salt concentration. In particular, we showed 
that the magnitude of the electrocaloric effect should decrease with increase in the initial 
temperature and increase with an increase in the salt concentration. Due to the general 
nature of the field theory approach Q| to tackle curved interfaces, polymers, multivalent ions 
etc., our work opens up a new area of theoretical research focused on the rational design 
of electrocaloric fluids. Furthermore, we have shown that the held theory approach stays 
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robust for high surface potentials and the other approaches such as the PB and MPB are not 
reliable. This particular feature of the held theory is quite important for energy harvesting 
technologies based on electrochemical capacitors and supercapacitors. 
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SUPPORTING INFORMATION: THEORY 


We consider two parallel plates separated by distance L and immersed in an electrolyte 
solntion containing Ug solvent molecnles, positive and n_ negative ions. The plates are 
assnmed to have uniform snrface charge densities (nnmber of charges per nnit area), ai and 
(72 (in nnits of electronic charge, e) and the corresponding snrface potentials are Vo,i and 
Vo 2 , respectively. It is to be noted that snrface potentials and charge densities are related to 
eachothe.bye.ect..„stat.cbo„„daryc„„ditio„sQ and depend on the mechanisms by which 
the plates acqnire the snrface charge. These relations can be formally derived by considering 
different mechanisms for charging l(] . 

Molecnlar volnmes of the solvent, positive and negative ions are taken to be Vs,v+ and 
n_, respectively. We are interested in nnderstanding the effects of dipolar interactions and 
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finite ion sizes on the thermodynamics of double layer. For such purposes, we seek a minimal 
model that can capture the underlying physics. In this work, we have studied a minimal 
model based on treating each solvent molecule as an electric dipole of length ps occupying 
molecular volume Vg- Also, the positive and negative ions have molecular volumes of u_|_ and 
V-, respectively. Finite polarizability { 2 , 3| of ions and solvent molecules are not considered 
in this work. However, the current formalism can be extended to take into account the 
effects of polarizability. 


The canonical partition function for such a system is written 


Z = 


as 


I Lr. , 

j=±,s a=l 


a=l 




( 1 ) 


\j = ±,s 


where rj^a is the position vector for the a*™ particle of type j and Uq, is the unit vector 
quantifying orientation of solvent dipole. The Hamiltonian is written by taking into 
account the contributions coming from ion-ion, ion-dipole and dipole-dipole interactions. 
Short range interactions between ions and solvent molecules are ignored in the minimal 
model studied here. Pj(r) represents microscopic number density of the particles of type j 
at a certain location r dehned as 


^ = s,+,- 


( 2 ) 


0 = 1 


The Hamiltonian for the ions and dipoles can be written as[ 

Ibo 


Ills 


H = 


dr / dr — 


pe{r) - Vr.P(r) 

Pe(r')-Vr'.P(r') 

r 

- r'l 


(3) 


where Ibo = /ineoksT is the Bjerrum length in vaccum and /3e(i') is the charge density 
(in units of e), given by Pe(r) = ZjPj{r) -|- ai6{x — Xi) -|- a 2 S{x — X 2 ), Zj being the 

valency (with sign) of ions of type j and |x 2 — Xi| = L is the distance between the plates. 
Also, P{r) is polarization density of dipoles (in units of e) at location r, given by 


P(r) = pg'^6{r- To) Uq 


(4) 


0=1 


Field theory in the canonical ensemble 


A held theory for the system described above can be constructed following a standard 


protocol 


3 . We start from the electrostatics contributions to the partition function. For 



















13 


the electrostatics contribution to the mrtition function written in the form Eq. [3l we use 
Hubbard-Stratonovich transformation!^ so that 

1 


exp 


-H 




exp 


i / dr pe(r) - Vr.-P(r) ^ ^/>(r) + 


Stt/ 


Bo 


dr-?/>(r)VrV’(r)) 


(5) 


where iVp is a normalization factor, given by 


= D [ijj] exp 


dr■ip{^r)Vl^p{r)) 


( 6 ) 


SttIbo 

Using this transformation and writing the local constraints (represented by delta functions) 
in terms of functional integrals using 


^ ^ M exp -i drr]{r) < ^ Pj{r)vj - 1 

r \j=±,s J I li=±,s 


(7) 


we can write the partition function given by Eq. [T] as 

1 


Z = ^ J D[iIj] J D[t]] exp 


H 


ksT 


( 8 ) 


where 


H 


SttI 


Bo 


dr^jJ{r)'\/l^jJ{r) — i / dr? 7 (r) + ai / dr||i'0(r||,Xi) 


+CT2 / drlli^p{rll, X2) - ^ {nj InQj {^jJ,T]} - Innj\} 

3=±,s 


(9) 


and we have used the notation r = {x,y,z) = (x,r||) so that ry denotes in-plane vector 
parallel to the plates. Qj is the partition function for particles of type j, given by 


Qj=±{i’,v} = J drexp[-iZj3/j{r) - ivj7]{r)] 

Qs{i>,r]} = dr duexp[-ipsU.'\/ri>{r) - ivsr]{r)] 


( 10 ) 

( 11 ) 


In the following, we use the saddle-point approximation to estimate the functional inte¬ 
grals over 'ijj and r]. An equivalent calculation in the grand canonical ensemble is presented 
in the Appendix A. 

















Saddle-point approximation: self-consistent equations, free energy and chemical 

potentials 
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The saddle point approximation with respect to rj and 'll! gives two non-linear equations. 
At the saddle-points, both rj and V’ turn out to be purely imaginary. Writing iri{T) = 77 * (r) 
and i'ip{v) = ' 0 *(r) at the saddle point and definmg densities of ions and solvent molecules 
via 


Pj=±(r) = 

P^(r) = 


n 


Qj h*} 

An'Us 


^ exp [—Zj'il!*{r) — Vjri*{r)] 

sinh ps| Vr' 0 *(r) 


Qs {-0*, V*} 

the equations at the saddle point are given by 


exp [-Vsr]*{r)] 




t=± 

V, ■ [e(r)Vri/>*(r)l = -4iriBoPe(r) 


( 12 ) 

(13) 


(14) 

(15) 


SO that the local charge density (pe(r)) and dielectric function (e(r)) are given by 


Pe(r) 

e(r) 


ZjPjir) + ai6{x - xi) + a2S{x 

i=± 


1 + inlBoPlPsir) 


L [p^|VrV’*(r)|] 

Ps\Vr^*{r)\ 


X2) 


(16) 

(17) 


where L{x) = coth(a;) — 1/x is the Langevin function. Corresponding Helmholtz free energy 
(F) is given by the approximation F/ksT = — InZ ~ H*/kBT = F*/kBT so that (cf. Eq. 

E) 


^ “ y drp*{r) + ai j dr|| 7 />*(r||, iri) 

+^2 f dry'll!*{rii, X 2 ) - {nj\nQj{'il!\p*}-Inrijl} (18) 

j=±,s 

Eq. [18] can be rewritten after eliminating rij using Eqs. [TH and [151 Furthermore, using 
the Stirling approximation In n\ ~ n In n — n, Eq. [18] can be written as 


F* 


dYpe{r)lp*{Y) + 


1 


drV'*(r)V^7/>*(r 


- / (irp5(r)ln 


dvr 


Stt/bo j 
sinh ps|Vr' 0 *(r)| 
p,|Vr^/'*(r)| 


+ Y ^^F(r) [Inpj(r) - 1] (19) 

j = ±,s 
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For study of opposing double layer systems in equilibrium with an electrolyte solution, 
chemical potential is determined by conditions in the solution far from the plates. In order 
to £x the chemical potentials by specifying different conditions in the solution far from the 
plates, we rewrite the above equations in terms of chemical potenials. An approximation 
for the chemical potentials (fij) of different species can be derived from Eq. [18] using the 
thermodynamic relation pj = {dF/dnj)^ ~ {dF*/dnj)Q = being the total volume. 

Using the Stirling approximation Inn! ~ nlnn — n, the chemical potentials within the 
saddle-point approximation are given by 


^j = ±,s 

ksT 


= In 


Un 




Using Eq. [201 Eqs. [12] and [13] can be written as 


( 20 ) 


Pi=±(r) 


exp 



dvr exp 



z,rij) - 

- Vsri*{Y) 


Vjr]*{Y) 

sinhps|VrV^*(r)| 

Ps|VrV^*(r)| 


( 21 ) 

( 22 ) 


Chemical part of the free energy: charging the electrodes and adsorption-desorption 

electrochemical equilibrium 


The free energy (cf. Eq. [T^ for the two opposing double layer system is obtained for a 
given surface charge density of the plates and has the charged plates at given surface poten¬ 
tials (in vacuum) as the reference frame. This can be easily seen by putting pj=±^s = 0 in 
Eq. [Ulsothat F^/ksT {pj=±,s = 0} = (cti/2) / dr||V'*(r||, aii) -F (ct 2/2) / dr||-?/’*(r||, a: 2 )- This, 
in turn, means that Eq. [19] doesn’t include the work done (typically by an external sourc^ 
in charging the two plates at a separation distance of L = |xi — a; 2 |. This contribution [3, 
to the free energy is 

P r /‘(T1 r rao 

^ chem 


ksT 


dr|| J {a'} 


drw 


da''ilj*{Y\\,X 2 ) {a'] (23) 


Evaluation of the right hand side in Eq. [23] requires specification of the mechanisms by 
which the plates acquire their charge. In the following, we consider the specihc case when 
plates are kept at constant surface potentials. 
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One dimensional model: plates at constant surface potentials with symmetrical ions 

and solvent molecules 


If the densities far from the plates are known to be pj^b corresponding to spatially nniform 
'^*(r) = and 'r]*{r) = pp then Eqs. ED and [22] can be written as 


Pi=±(r) = Pyfecxp [-Zj {V^*(r) - ^pp} - Vj {r/*(r) - pp}] 

/X r r ^ *11 sinhp,|Vr^/’*(r)| 

Ps{r) = Ps,bew[-Vs{v (r)-%}] 


(24) 

(25) 


Ps\Vri>*ir)\ 

For two parallel plates, saddle point valne of pj varies only along the direction perpendic- 
nlar to the charged snrface (taken to be along x-axis) so that = ijj*{x), p*{r) = p*{x). 

Fnrthermore, considering the case of symmetric ions and solvent molecnles so that Vj=±^s = 
and = jZcj so that Pj=±,b = Pc,b, we can eliminate p* nsing Eqs. [HI ED and 


and write Eq. [15] as 


A 

dx 


e(x) 


d'il)*{x) 

dx 


= -iTllsoPeix) 


(26) 


where the local charge density {pe{x)) and dielectric fnnction {e{x)) are given by 
Pe{x) = \Zc\[p+{x) - p_{x)] + aiS{x - xi) + a 2 S{x - X 2 ) 

/ x ^ Pc,bew[-\Zc\{P’*{x) -V’fe}] 

/X _ Pc,bew[\Zc\{'ip*{x)-ipl}] 



Ps,b 

sinh ps 

dp*{x) 

dx 

L 

Ps 

1 dp*(x) 

1 dx 



[r{x) - 

■ Ps 

dp*{x) 1 

dx 1 

Ps\ 

d'il)*{x) 1 

dx 1 


(27) 

(28) 

(29) 

(30) 


SO 


that 


sinhn 

^cosh [\ZP\ {^p*{x) - 'ipl}] (31) 

Ps\—Q^\ 

and [psfi + 2pc,b] = 1. It is to be noted that solvent density is given by 


fU = p.xa' 

'■ ’ dx 


sinhp,|^^| 


and 


. _ _ p^ _ smhp^p 

~ f N dib*(x\^ „ idp*(x) 

satishes the incompressibility constraint g Pj{x)a^ = 1. 


( 32 ) 
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Free energy within saddle-point approximation : adiabatic changes 


Changes in entropy {AS) can be readily calculated from the corresponding free energy 
changes (AF) and the thermodynamic relation AS = — Free energy of the double 

layer system {F^i) is the sum of electrostatic contributions approximated by F* and the 
chemical part given by Superscript -k implies the use of saddle-point approximation 

(mean-held like treatment) in estimating the free energy. In particular, assuming lateral 
homogeneity, for plates (at known surface potentials) separated by distance L having surface 
area A each, F* and are given by 

F* . 1 


AkeT 


/ dxpe{x)i)*{x) + 

/o otiLbo Jo 

.L r sinhp, 


dx'ip {x) 


dxps{x) In 


'0 


dvr- 


dx 


I dp*{x) I 
Ps\ \ 


dx"^ 

J=±>« 


and 


F 


chem 


= -ai'il:*{xi) - a2'ip*{x2) 


dxpj{x) [lnpj(a;) — 1] (33) 


(34) 


AksT 

In order to compute the electrocaloric effect, free energy changes with respect to the system 
in the absence of applied electric held are desirable. In the absence of applied electric held 
(i.e., when ai = ct 2 = 0 and considered as the reference state), the same number of ions and 
solvent molecules are homogeneously distributed in volume = AL so that free energy of 
the reference state becomes 


ALksT 


chem 


ALksT 


= 2pcfi [In Pc,b - 1] + Ps,b [hi Ps,b - 1 - In dvr] (35) 


O'! = 0 - 2=0 


where, we have used the constraint A dxpj^^ = pj^^ for equating the number of ions and 
solvent molecules in the absence and presence of applied electric held. Using these equations, 
the free energy changes (AF*) due to the application of an electric held can be written as 


AF* 


^ dl 


AksT AkuT 




~ “ rr2ili'{x2) 


(36) 


Spatially uniform dielectric : Poisson-Boltzmann (PB) and modified 
Poisson-Boltzmann (MPB) approaches 


In the limits of small surface potentials so that '4)*{x)—'ipl —>■ 0 and weak coupling limt for 
dipoles, dehned by Ps | I 0 , the dielectric function given by Eq. 130) becomes spatially 
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uniform so that 


47r 

€{x) = €h ^ 1 + -^IboPIPsX 


(37) 


Physically, this means that solvent density is spatially uniform in the limits of small surface 
potentials and weak coupling limit for dipoles so that Ps{x) = ps,b as evident from Eq. 
32l It is to be noted that the quantity / is taken to be unity in these limits and leads 
to the standard Poisson-Boltzmann results pioneered by Verwey and Overbeek[^. Another 
somewhat recent development (so called modified Poisson-Boltzmann (MPB) approach Qj) 
is to consider the case of uniform dielectric but include steric effects in the calculations of 
charge density by taking 


/ 1 


d'ip*{x) \ _ 


dx 


= - ifb} = 1 - ao + ttocosh [\Zc\ {'f*{x) - V’h}] 

(38) 


, where cxq = 2pcfiaf‘ is the packing fraction of ions in the bulk. Although it seems incon¬ 
sistent to ignore and retain functional dependence of a particular quantity such as / while 
considering different physical quantities such as dielectric function and charge density, the 
MPB approach has been quite successful in predicting qualitative features of the double 
layer capacitance. Nevertheless, the MPB approach leads to semi-analytical predictions for 
the electrostatic potential and the free energy, as described below. 

With the approximations described above, Eq. [26] can be readily integrated over x (after 
multiplying by on both sides). In particular, we obtain a self-consistent equation for 

dip*(x) 

dx 


2 dx ^ ^ -A}-M (39) 

where A is an integration constant, which is determined below and the effects of surface 
charge densities (cti,(T 2 ) appear in the form of boundary conditions. Using Eq. 15^ and 
equations at the saddle-point, the free energy changes of the double layer system, defined 
by Eq. [36l can be written as 


AkeT AkeT ^ ^ ' ' 

where F^i^mpb = Fmpb + ^chem ^^^1 is the approximation for Eq. [33] obtained using 

Eq. [39] and is given by Eq. [M] We must point out that in obtaining Eq. [lO] we 
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have retained functional dependence of the solvent density on /mpb through Eq. W2\ and 
used the incompressibility constraint. 

In the following, we consider two cases of non-overlapping and overlapping double layers 
and eliminate A from Eq. 001 In the case of non-overlapping double layers, 'ip*{x) becomes a 
non-monotonic function of x with a minimum at a; = Xmin- Integrating Eq. [39] over x with 
the limits Xi and X 2 , we obtain js] 




+ 


dip* 


[Ih flMPB (x) “ Ip^} “ A] J [Ih fjMPB (x) “ Ip^} “ A] 


Stt/ 


Bo 


-,1/2 




(41) 


Similarly, for the case of overlapping double layers so that ip*{xi) > ip*{x 2 ), we obtain 


cip*{xi) 


dip* 


/p*(x 2 ) [In Jmpb {iJ*{x) - ipp} - A]^/^ 


Snl 


Bo 




1/2 


(42) 


Eqs. mi and 02] allows us to eliminate A from Eq. 00] and write it as 


A F* 
mPB 


eh 


-g{ip*{xi),ip*{x 2 )} = 


‘^Pc,b\Zc 


AksT Y 2'kIbooP'’ ^ ^ 

where we have defined = S'kIbo\Zc\^ pc,blen- Also, 

rP*{^k) j - 

g{ip*{xi),ip*{x2)} = ^ / dlp^J\TlfMPB{'ip*{x) -ipl) 


■g{ip*{xi),ip*{x 2 )} (43) 


(44) 


for the non-overlapping double layers and 


g{ip*{xi),ip*{x2)} 


l<p*{xi) 

Ill)*{x2) 


diipXn Smpb {'ip*{x) - V'P 


(45) 


in the case of overlapping double layers. 

Changes in entropy {ASI^pb) can be readily calculated using Eq. 03] and the thermody¬ 
namic relation AS = — so that 

^^MPB _ _ ( ^^^MPB ^ _ 6Pc,fcI Zc I 

AkB V AkBdT y/2B?a~o^ 

where we have dropped explicit functional dependencies of g on ip*{x) — ipp for convenience 
in writing. It is interesting to consider the limit of dilute solutions so that ckq —)■ 0 and this 
limit is the same as the standard PB approach. In this limit, for non-overlapping double 


1 + 


3gdT 


(46) 
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layers, 'ip*{xmin) = V’b and A = 1 (due to the fact that d'ip*{x)/dx = 0 ai x = Xmin in Eq. 
1391) . This leads to 


AkeT 


4p, 


c,6 


K 


^ [2 - 2cosh(\Z4 


(47) 


k=l,2 


i.e., the total free energy change is the sum of changes in the individual double layers Sj. 
This leads to entropic changes given by 


Aks 


^Pc,b 


K 


E 


fc=l ,2 

1 — sech 


3 — 2 cosh 

Z. 




- V’H ) - sech 

T den 


{'lp*{Xk) - V’b} 


ehdT 


(48) 


NUMERICAL METHODS 


We have solved the set of equations (Eqs. [26]- 13^ numerically after rewriting Eq. [26] in 
the form 

di>*{x) _ d‘^ii*{x) 1 de{x) di/j*{x) AtiIbo , . . . 

dt dx‘^ e{x) dx dx e{x) ^ 

where f is a hctitious time. A steady state solution of Eq. 09] is obtained by using the 
extrapolated gear scheme and using size of ions a to obtain dimensionless length variables. 
Time step of 0.0001 is used to integrate Eq. 09] with L/a = 20 — 40 (depending on the 
value of pcfi) and 1024 grid points. Convergence of the numerical solution is checked by 
computing free energy changes between two consecutive time steps and the changes less 
than 0.0001 are used to set the tolerance criteria. These equations are solved for non¬ 
overlapping double layer systems so that one of the surfaces has the known surface potential 
while the other is grounded (i.e., surface potential is zero). The temperature is changed by 
varying Ibo and the free energy changes (in units of AksT/a?) are computed using Eqs. [33l 
1341 [35] and [36] In computing the electrocaloric effect, we have made use of the fact that 
the held variable 'ip*{x) in the theory is the electrostatic potential (in units of ksT/e) at 
location x. For example, '^*(0) = eVby/fcsT for the single double layer system studied in 
this work. Numerical estimates for the surface charge densities we obtained by the relation 
cr = -[e{x)/A'KlBo){d'4^*{x)/dx\^=Q. 


















21 




(a) (b) 

FIG. 4: (a) Effects of the bulk salt concentration on the free energy changes {AF = AF*) of the 
double layer at To = 303 K. (b) Comparisons of electrostatic potential profiles from the 

MPB approach and numerical calculations at Cs = 0.1 M. 

RESULTS: ANATOMY OF THE DOUBLE LAYER 

Anatomy of the double layer is determined by the electrostatic potential prohle. As 
the comparisons between the PB and MPB approaches are well-known j^, we only show 
comparisons between the MPB and our numerical calculations in Figure ID^b) for low and 
high surface potentials. It is found that the MPB and numerical results are in excellent 
agreement at eVo/^s^o = 1 showing exponential decay appearing as linear on semi-log 
plot, as expected. In contrast, the electrostatic potential prohles differ near the surface (for 
x/a < 2) at eVo/ksTo = 10, which are responsible for differences in free energies predicted 
using the MPB approach and the numerical calculations (cf. Figure 3(a) in the main text). 

The differences in the electrostatic potential near the surface show up in plots for surface 
charge density (ai = cr) as a function of applied surface potential (Figure 3(b) in the main 
text). The structural changes resulting from an increase in the surface potentials are shown in 
Figure [5l In particular, an increase in surface potentials leads to an increase in the volume 
fraction of counterions (anions in this case) near the surface at the expense of excluding 
coions and solvent molecules. However, further increase in the surface potential (e.g., see 
plot for eVo/ksTo = 9 in Figure [5](c)) leads to increase in solvent volume fraction near the 
surface at the expense of exclusion of counterions and coions. This is expected from the 
expression for solvent volume fraction, Eq. [32], leading to higher volume fraction of solvent 





22 



(a) 




(b) 



x/a 


(c) (d) 

FIG. 5; (a) The dielectric function, (b) electrostatic potential, (c) solvent and counterion (anion) 
densities, and (d) co-ion (cation) densities at different surface potentials are shown for bulk salt 
concentration of 0.1 M at Tq = 303 K. 


in regions having strong electric fields. Also, such an enrichment o: 


so 

id, 


vent in regions of 


3 . 


Furthermore, 


strong electric helds is in agreement with previous theoretical works 
the electric held dependent sorption of water on the AFM tips has been used to modulate 
friction at the nanoscale 12]. 
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APPENDIX A : Field theory for double layer systems in the grand canonical 

ensemble 

For study of a double layer, grand canonical partition function can be constructed and is 
given by r = s 


so that 



D [rj] exp 


Hg{'lp,V} 

ksT 


(A-1) 


j dv%lj{jc)Vl'il){Y) - i j drr7(r) + cri j dr\\i'ip{r\\,xi) 

+^2 f driii'ip{r\\,X2) - ^ Q^{^,r]) (A-2) 

i=±.s 

where we have used Eqs. |8]and |9]for the partition function in the canonical ensemble. 

The saddle point approximation with respect to rj and V’ gives two non-linear equations. 
At the saddle-points, both r] and ifj turn out to be purely imaginary. Writing ir]{r) = ri*{r) 
and at the saddle point, the two equations are given by 


'^VjPj{r)+Vsps{T) = 1 (A-3) 

i=± 

Vr ■ [e(r)VrV'*(r)] = -^nlsoPeir) (A-4) 


where we have defined 

pj=±{r) = exp 


k'j 

ksT 


Zji^*{r) - Vjp*{r) 


Ps{r) = dvr exp 


Ps 


ksT 


- VsP*{r) 


sinh ps|Vr'0*(r)| 
Ps\Vri’*{r)\ 


(A-5) 

(A-6) 


so that Pe(i') = J2j=± ZjPj{r) -|- aiS{x — xi) -|- ct2 (5 (x — 2:2) and the local dielectric function 
is given by 


e(r) = 1 +A-nlBoPsPs{r) 


2 / xA [p.|Vr^/’*(r) 


(A-7) 


p^|Vr^/’*(r)| 

where L{x) = coth(a;) — 1/x is the Langevin function. Corresponding approximation for the 
Gibbs free energy is given by 

1 


ksT 


= ai dr||'0*(r||,a;i)-h 0-2 / dr 11-0*(ry,2:2) + 


SttI 




-J dr?7*(r) - ^ j drpj(r 


Bo 


j=±,s ' 


(A-8) 

















Using Eqs. IA-31IA-41 IA-51 and IA-61 it can be shown that H* and F* given by Eq. [19] are 


related by 



(A-9) 


in accordance with the thermodynamic relation that the Helmholtz free energy is the Gibbs 
free energy plus chemical potential times the number of particles. 
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